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We study KPZ surfaces on Euclidean lattices and directed polymers on hierarchical lattices subject 
to different distributions of disorder, showing that universality holds, at odds with recent results 
on Euclidean lattices. Moreover, we find the presence of a slow (power-law) crossover toward the 
universal values of the exponents and verify that the exponent governing such crossover is universal 
too. In the limit of a 1 + e dimensional system we obtain both numerically and analytically that the 
crossover exponent is 1/2. 
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The problem of Directed Polymers in Random Media 
(DPRM) (!]] has attracted much attention in the last ten 
years, both as a paradigm in the area of disordered sys- 
tems and for the richness of its connections with other 
systems, in particular noisy surface growth governed by 
the Kardar-Parisi-Zhang (KPZ) equation jg||. Both 
problems show space and time scaling behavior, and the 
connection between the two manifests through a corre- 
spondence between their exponents. Within the KPZ 
context, it has been recently proposed that the model 
could be non-universal Q] : The exponents characterizing 
the surface growth depend on the details of the driving 
noise, at least for a substrate dimension d > 2. These re- 
sults have been obtained via a lattice formulation of the 
KPZ equation where the strong coupling limit (which 
is the non trivial regime of surface growth) is shown to 
be completely equivalent to the ground state problem of 
DPRM's. Thus, we can expect non-universality in the 
DPRM context as well. It is therefore interesting to re- 
examine this issue both within the KPZ context and for 
DPRM's. 

Before proceeding further, it is useful to recall the defi- 
nition of the exponents within the two contexts, and their 
relations. Starting from a flat substrate of characteristic 
linear size L, the roughness of a KPZ surface grows ini- 
tially as W(t, L) ~ t@; at longer times it saturates, and it 
scales with L as W(t, L) ~ L x . The characteristic time r 
between the two regimes scales with the size of the system 
as t ^ L z . These exponents are not independent: the re- 
lation (rooted in the Galilean invariance of the KPZ equa- 
tion) x + z = 2 holds in every dimension |^,|3| . Moreover, 
consistency imposes x = Z P- Therefore there is just a sin- 
gle independent exponent. In the language of DPRM's, 
the exponent /3 governs the fluctuations of the ground 
state energy, AE = y/< (E GS - < E G s >) 2 > ~ t u with 
u> = ft (here we use the greek alphabet letters commonly 
used in the literature); the transverse wandering fluctu- 
ations of the ground state polymer are governed by the 
exponent ( = 1/z, Al ~ tS. The relation uj = 2£ — 1 
holds in every dimension (it is related to x + z = 2). 

The exponents are known exactly only for d = 1: 
P = 1/3 and z = 3/2 (w = 1/3, C = 2/3). At present 



there are no exact solutions for d > 2, and our knowledge 
of the exponents relies on numerical simulations. Yet, the 
situation is not completely clear even numerically: As it 
has been pointed out recently , in the context of surface 
growth (KPZ equation) different distributions of disorder 
seem to give different values of the exponents. This re- 
sult has been interpreted as a case of non-universality. 
Indeed, using a distribution of the energies as 

pW = ^(i-Hn ze(-l,l), a>-l , (1) 

in pj it was found that the values of the exponents de- 
pended on the value of a in ([[J). In the exponents 
were calculated only for d > 2, assuming that universal- 
ity holds in d = 1, where the exponents are exactly known 
(it is worth mentioning that the exact knowledge of the 
exponents is based on taking a Gaussian distribution of 
disorder). 

To settle this problem, we have performed simulations 
of surface growth on d = 1 Euclidean lattices. We mea- 
sure the (3 exponents ruling the growth in time of the 
roughness W(t) starting from a flat substrate (the other 
exponents can be obtained from the above mentioned re- 
lations). We used three different values of a, namely 
a = 0.5 ,0.75 ,0.9 (in each case we also verified that 
we are still far from saturation). As it can be seen from 
Fig.|l|, a power law with exponent 1/3 (the theoretical 
one) is not suited to fit the numerical power laws ob- 
tained. A naif fit would give instead = 0.28 , 0.26 , 0.24 
respectively. Even fitting an exponent on the last four 
points with a = 0.5 would give (3 = 0.30 (results for the 
most commonly used case, the uniform distribution with 
a = 0, are not given since the fitted exponent is /3 = 0.31 
and the crossover is less evident). To understand whether 
this is a case of non-universality or of crossover, we also 
analyze the running exponents (see Fig.^Ja). If the numer- 
ical points shown in Fig.[l] are taken at fixed time ratios 
tn+x/tn = then the running exponents are defined as 



Pn = log, 



W(t n+1 ) 
W(t n ) 



(2) 



where, in the present case, r = 3/2. We find that the 
running exponents approach their universal value j3 = 



1 



1/3. In a log- log plot (Fig J^b) it is then easy to see that 
such an approach is ruled by a power-law with a universal 
(independent of a) exponent 7 = 0.23 ± 0.02. 

The results obtained are amenable of different interpre- 
tations: either the d = 1 case is non universal too, against 
naive expectations, or it is indeed affected by extremely 
slow crossovers. The analysis of the running exponents 
seems to sustain the latter hypothesis. The behavior in 
higher dimensions is unfortunately much more difficult to 
extract. Indeed, we have also performed simulations in 
d = 2, but the results, although similar, are somehow less 
conclusive due to numerical limitations: in order to get 
reliable statistics in d = 1, each numerical point is aver- 
aged over 10000 independent realizations. Moreover, uni- 
versality in the crossover is evident only for times > 100. 
Matching both conditions for d > 2 goes beyond the ca- 
pabilities of our present computing facilities. 

To investigate the high dimensional behavior of this 
crossover we resort then to analytical and numerical cal- 
culations on hierarchical lattices, where the same behav- 
ior emerges. 

The main idea of hierarchical lattices is that it is possi- 
ble to build them iteratively, given a fundamental bond- 
block transformation, where a bond is substituted by a 
block as in Fig.|. The inverse process can be seen as a 
coarse graining transformation. Indeed, using this trans- 
formation, it can be shown that real-space renormaliza- 
tion becomes exact on hierarchical lattices. We exploit 
this property to write the corresponding rcnormalization 
for the ground state of polymers on hierarchical lattices 
with bond disorder taken from a given distribution. 

Given the ground state energy distribution at a cer- 
tain renormalization step n, it is possible to compute the 
distribution of the ground state energy at step n + 1 via 
the equation 



uj n = log 2 



Cn+l 



Or, 



(4) 



P n+ i(x) = bQ n {x) 



Q n {x')dx' 



b-1 



(3) 



where b is the number of sides as from Fig.|| Q n (x) 
is the convolution of two probability distributions P n (x), 
Qn(x) = P n (x)*P n (x). The r.h.s of Eq.(||) represents the 
probability that one of the b sides has an energy x, while 
all the others have an energy greater than x. Due to the 
choice of the minimum energy, this is also the probability 
that the ground state energy is x at step n + 1. 

The number of sides b can be related to a fractal di- 
mension of the lattice. Indeed, at every renormalization 
step we rescale the length of the lattice of a factor 2, and 
the volume of a factor 2b. Therefore the dimension of the 
system can be related to b via the formula D = 1 + log 2 b 
(or, correspondingly, d = D — 1 = log 2 b). 

Given Eq. (||) , the calculation of the exponent u reduces 
in principle to a numerical iteration of (|3|) given a starting 
distribution Pq(x). The ground state energy fluctuation 
exponent can be computed as 



where a n is the variance of P n (x). Asymptotically, the 
running exponent u) n tends to a constant lOqo. Possi- 
ble non-universalities should then emerge using differ- 
ent distributions. Yet, calculations of show that 
different distributions give the same asymptotic values 
up to b = 20 (corresponding to a substrate dimension 
d = 4.32...). 

As a byproduct, we see that hierarchical lattices of 
fractal dimension D = d + 1 give values of the to expo- 
nents close to the results on Euclidean lattices of the same 
dimension (D = 2, uj hier = 0.30(1), w euc ; = 1/3; D = 3, 
Wfcier = 0.22(1), uw = 0.24(1); D = 4, u) hier = 0.15(1), 
uj euc i = 0.16(1)). 

Although universality holds for directed polymers on 
hierarchical lattices, we find that the running exponents 
uj n show a power-law crossover toward their asymptotic 
value Woo (see Fig.^Ja): 



-7(n+l) 



(5) 



2«+i being the length of the lattice. 

The power-law exponent 7 does not depend on the 
details of Pq(x), that instead are responsible of the am- 
plitude of the crossover. We used Gaussian distributions 
and distributions of the form ([!]) , as considered in Q . 

In Fig.|| we show the running exponent approach to 
Woo for different values of b (all the numerics in this case 
have been done using Gaussian distributions; as already 
explained above, universality of the crossover exponent 
is verified). We find that, as b — > 1 then 7 — > 1/2. 

We verify that indeed 7(6 = 1 + ) = 1/2 via a b = 1 + e 
expansion of Eq. (|J) on the same lines in Ref. Q : there a 
Gaussian microscopic distribution Pq(x) was used, and 
no crossover effects were discovered. We use instead 
gamma distributions of the kind P(x; v) = Yjy ) x w ~ 1 e~ x , 
defined for x > 0: these distribution are partially stable 
under convolution; indeed, P(x; v) * P(x; v) = P(x; v' = 
2u), that is, the convolution of two gamma distributions 
is still a gamma distribution (although characterized by 
a different exponent v'). 

We assume that the distribution at iteration n is given 
by 



P n (x)=P°{x) [l + e$„(z)] 



(6) 



where P®(x) is the distribution at iteration n for the 1- 
dimensional system (b = 1), simply given by Pn(x) = 
P^i{x) * P°_i(a;) = The perturbation $„(x) 

must satisfy the relation 



$ n (x)P°(x)dx = 



(J) 



We expand Eq.(|3|) for small e, keeping only terms of 
order e: 



2 



Q° n (x)+ e Q° n (x)$ n+1 (x) = Q°Jx) + 

+ e {Q°(z)+Q°(x)ln / g°(x')^' + 

+ 2P%(x)*[P°(x)$ n (x)]}. (8) 

The linear terms in e are a recursion relation for the 
perturbation &(x), 



Q° n (x)$ n+1 (x) = Q° n (x) + Q° n (x)ln / Q° n (x')dx' + 

J X 

+ 2PZ(x)*[P°(x)$ n (x)]. (9) 

All the terms in (||) are proportional to Q° (x) except for 
the last term on the r.h.s. We deal with this last term 
with the ansatz equation 

P°(x) * [P°{x)4^(x)] = \ n , s Q° n (x)<Pn+lA X ) ( 10 ) 

whose purpose is to extract from the l.h.s. expression a 
term proportional to Q„(x) that can be then simplified 
in (||) . If among the solutions of Eq. ( |l~o| ) it is possible to 
find sets {<fi n ,s} that are complete and orthonormal for 
any n, then we can write & n (x) = J2 s a n,s4'n,s(x) and 
@ becomes a recursion relation for the a n coefficients. 
Any function of the form 



(-i) s rK) 



! I/ T(v n + s) P°{x;v n ) dx s 



satisfies Eq.(lO) with "eigenvalues" 



' T(v n + s)T{u n+ i) 
T{v n )T(2v n + s) ' 



(11) 



(12) 



Eq. (|1 1[) is nothing else than Rodrigues' formula for the 
generalized Laguerre polynomials, properly normalized, 
that are known to be orthogonal and complete with 
weight Pn{x) 0. The first three such polynomials are 

<t>n,o(x) = 1 , 

4>n,l( x ) = (x~ V n )l^fv^ , 

<$>n,i(x) = {x 2 - 2(z/„ + l)x + v n {v n + 1))/ \j2v n {v n + 1) . 

(13) 

Substituting the series expansion of < &„(x) in (^), it is 
possible to read the iteration equations for the coefficients 



a n + l,0 — 1 + Knfi + 2A ra ,oOn,0 

a-n+i.s = Kn,s + 2X niS a niS s > 



(14) 



with 



Kn,. = J dxP° +1 (x)(fr n+ i tS ln I P° +1 (x')dx' (15) 



It is straightforward to show that Kq = — 1; moreover 
from Eq. (Q) we find ao,o = an d therefore a„ t o = for 
any n, as required. It is also possible to show that 2\ nyS < 
1 if s > 2 (indeed the leading behavior for large v n , that 
is for large n, is 2A„, S - 2 1 ~ s / 2 [l + (s - l/2)s/2v n ] 1 / 2 ). 
Therefore all the iterations equations (|l4|) with s > 2 
converge and are asymptotically irrelevant for the scaling 
behavior. We are left with s = 1,2. Indeed a n ,i and 
fln, 2 are the only coefficients that are relevant for the 
computation of u) n from (|J). 

We can then compute the variance a\ as 



al = < x 2 >„ - < x >l = 



■ n + e{^/2v n (y n + l)a„ >2 + 2y/vZa n ,i) (16) 



= v. 

We can write then 

T 2 



(T 



rt+1 
0-2 



1 



V2^„ 



(17) 



where we keep only the leading and next-to-thc-lcading 
terms. From ( |l7j ) the order e correction to the u> exponent 
can be computed. Indeed, (J 2 + i/<7 2 = 2 2aj ™, and since 
u n+ i = 2 n+1 VQ, we can write 



e hT2 K 



n,2 



(18) 



Furthermore, both K n> x and K n ^ have a power-law con- 
vergence to their asymptotic value, with exponent 1/2 
(this result can be obtained using a Gaussian approxi- 
mation in their evaluation). As a result the r.h.s. of 
converges with exponent 7 = 1/2 to Wqo. 

In conclusion, we have found numerically that both the 
KPZ equation on d — 1 Euclidean lattices and DPRM 
on hierarchical lattices show universal behavior, being 
the exponent ui independent on the details of the disor- 
der distribution. Yet, the presence of non trivial power- 
law (therefore " slow" ) crossover effects has been unveiled. 
Moreover, we have shown that the crossover exponent 7 
is universal. We have also shown numerically and ana- 
lytically that 7 — * 1/2 when b — ► 1. 

Since wc find the same qualitative slow crossover be- 
havior both on d = 1 Euclidean lattices and on hierarchi- 
cal lattices of any dimension, it is reasonable to think that 
the same behavior is present also for d > 2 Euclidean lat- 
tices, providing a way out of the non-universality claimed 
in|. 

We thank C. Tebaldi for useful discussions. 
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FIG. 3. Iterative procedure to build a hierarchical lattice: 
the bond on the left schematically represents the lattice at 
step n; then it is used to build the lattice at step n + 1, 
represented on the right. In the picture only three of the b 
sides are explicitly drawn. 




FIG. 1. Roughness of KPZ surfaces for exponents 
a — 0.5,0.75,0.9 in Eq.(l). The dashed line represents a 
power-law with exponent f3 — 1/3, the dotted lines represent 
temptative fits. Each numerical point is the average from 
10000 independent disorder realizations. 
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FIG. 4. a) Running exponents for 6 = 8 with different dis- 
order distributions, namely Gaussian and as from Eq.(4); b) 
Log-log plot of the difference of the same running exponents 
from their asymptotic value, clearly showing the power-law 
behavior with 7 = 0.80(2). The value b = 8 has been chosen 
for numerical stability. 
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FIG. 2. (a) Running exponents for the roughness as in 
Fig.l. (b) Log- log plot of the difference of the running ex- 
ponents from the theoretical value 1/3. The dashed line rep- 
resents a power-law with exponent 0.23. 
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FIG. 5. Log-log plot of the difference of the running ex- 
ponents from their asymptotic value for different values of b. 
The dashed line represents a power-law with exponent —1/2. 
All the data have been obtained using Gaussian distributions. 
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